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Molecular dynamics simulations were conducted to investigate the structural properties of melts of noncon- 
catenated ring polymers and compared to melts of linear polymers. The longest rings were composed of 
N = 1600 monomers per chain which corresponds to roughly 57 entanglement lengths for comparable linear 
polymers. For the rings, the radius of gyration squared, (-R 2 ), was found to scale as TV 4 / 5 for an intermediate 
regime and TV 2 / 3 for the larger rings indicating an overall conformation of a crumpled globule. However, 
almost all beads of the rings are "surface beads" interacting with beads of other rings, a result also in agree- 
ment with a primitive path analysis performed in the next paper 1 . Details of the internal conformational 
properties of the ring and linear polymers as well as their packing are analyzed and compared to current 
theoretical models. [DOI: 10.1063/1.3587137] 



I. INTRODUCTION 

Understanding the static and dynamic properties of 
ring polymer melts is one of the remaining challenges in 
polymer science. Unlike linear polymers, the topolog- 
ical constraints for ring or cyclic polymers are perma- 
nent and this affects both their static and dynamic prop- 
erties. For linear polymers the topological constraints 
imposed by the non-crossability of the chains (entangle- 
ments) force each chain to diffuse along its own coarse- 
grained backbone, the so-called primitive path, and this 
is well described by the reptation model of de Gennes 
and Edwards. For branched systems strands have to fold 
back in order to find a new conformation without cross- 
ing any other chain, resulting in an exponential growth 
of the longest relaxation time due to the entropy bar- 
rier of O (strand length) between different states^H A 
number of simulations and experimental results confirm 
this conceptP^. For both linear and branched polymer 
melts, it is the free chain ends that make the known re- 
laxation mechanisms possible. However, in the case of 
rings there are no free ends. This raises a number of 
unanswered questions regarding the motion and stress 
relaxation of ring polymer melts based on their different 
conformational properties, which will be discussed in the 
subsequent paper^. 

Understanding a melt of nonconcatenated and unknot- 
ted rings is also closely related to the problem of the 
compact globular state of an unknotted loop. Grosberg 
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et al. 14 were the first to hypothesize that the equilib- 
rium state of a compact unknotted loop is the so-called 
crumpled globule, in which each subchain is condensed 
in itself, and, therefore, the polymer backbone is self- 
similar with a fractal dimension of 3. The idea behind 
this prediction is that two pieces of the chain, each of 
them crumpled, act pretty much like two nonconcate- 
nated rings, and the latter system obviously experiences 
at least some volume exclusion for topological reasons (on 
this subject see recent work— and references therein). A 
direct test of this prediction was first made by Lua et 
alP^ by modeling the closed Hamiltonian loops on the 
compact domains of a cubic lattice. Some evidence of 
segregation between globules on the chain was observed, 
but overall the results were inconclusive because the sim- 
ulated chains were too short or the statistics were too 
poor. 

Squeezing topologically constrained rings against each 
other in a melt, or squeezing one unknotted ring against 
itself in a restricted volume, is also a problem of great 
potential significance in application to the DNA organi- 
zation in the cell nucleus. Chromatin fibers are packed 
in vivo at a rather high density, with volume fractions 
not dissimilar to those in a polymer melt. However, un- 
like a melt of linear chains, the different chromosomes in 
the nucleus do not intermix, but stably occupy different 
distinct "territories": the image in which each chromo- 
some is stained with a different color resembles a po - 
litical map of some continent 17 (see also reviews^^) . 
This makes an obvious analogy with the melt of non- 
concatenated rings; indeed, if the rings in the simulated 
melt are shown in different colors, the image of a political 
map emerged. This strongly suggests topological prop- 
erties of chromatin chains as a likely mechanism behind 
chromosome segregation, similar to segregation of non- 
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concatenated rings in the merP^n] Why the topology 
of the chromatin fibers is restricted is a somewhat open 
question, but it might be that the DNA ends are attached 
to the nuclear envelope, or simply that the cell lifetime 
is not long enough for reptation to develop^HSD Territo- 
rial segregation of chromosomes appears to be a common 
feature for the cells of higher eukaryotes, including hu- 
mans, but it is less well pronounc ed or ab sent altogether 
in lower eukaryotes such as yeastP^ZEHl. This observa- 
tion is at least qualitatively consistent with the fact that 
topological segregation of nonconcatenated loops is fully 
developed only when polymers are very long. 

In recent years, more detailed information on chro- 
matin fiber organization has become available, due to 
the advent of novel experimental techniques, such as Flu- 
orescence In Situ Hybridization FISH 29 and HiC, which 
is a generalization of 3C (Comprehensive Chromosome 
Capture^Pand 4C (the same on a Chip) 31 . These meth- 
ods allow one to probe the large scale features of the 
chromatin fiber fold in a single chromosome. In partic- 
ular, the HiC study of human genome folding 32 revealed 
the signature of a fractal folding pattern which was ex- 
plicitly associated with the crumpled globule organiza- 
tion of DNA originally predicted on purely theoretical 
grounds^. Specifically, the measurement indicated that 
the loop factor P(s) (the probability that two loci of ge- 
netic distance s base-pairs apart will be found spatially 
next to each other in chromosome) scales as P(s) ~ s~ 7 
over the interval of s, roughly, 0.5 Mbp < s < 7 Mbp, 
where the power 7 is 1 or slightly higher. This scaling ap- 
pears to be consistent with the crumpled globule model 34 . 
There are also many more indications of self-similarity, 
or scale-invar iance, in the chromatin folding (see recent 
review^ and references therein). Current biological the- 
ories about these issues also actively involve dynamical 
aspects, how this presumably crumpled self-similar struc- 
ture can move to perform its functions. To this end, 
there are some observations regarding the structural dif- 
ference between the active part of the chromatin fiber, 
currently transcribed, and the more densely packed, non- 
transcribing heterochromatirP 2 -'. Without going into fur- 
ther detail, suffice it to say for the purposes of the present 
paper that chromosome studies necessitate the further 
and deeper understanding of the melt of nonconcatenated 
rings, both in equilibrium and in dynamics, as a basic el- 
ementary model. 

Recent theories for melts of nonconcatenated rings or 
the globule of a single unknotted ring have been dom- 
inated by the idea that the topological constraints im- 
posed on a given ring by the surrounding ones, as well 
as the constraints imposed by the different parts of the 
same ring on each other, can be treated using a sort 
of a "topological mean field", namely, by considering 
a single ring immersed in a lattice of immobile topo- 
logical obstacles^MlS w ith none of the obstacles pierc- 
ing through the ring (such that the ring is nonconcate- 
nated with the whole lattice). The conformation of a 
ring in such a lattice represents a lattice animal with 



every bond traversed by the polymer twice (in opposite 
directions). We emphasize that the structure of branches 
in the lattice animal is completely annealed. In the off- 
lattice case, a similar annealed randomly-branched struc- 
ture is expected, in which the chain follows itself twice 
along every branch only down to a certain scale <i, sim- 
ilar to the distance between obstacles or the tube di- 
ameter in reptation theory. Such structure is reason- 
able as long as Ik <C d <C (R g ), where Ik is the Kuhn 
length and (R g ) is the average gyration radius of the loop. 
Based on the lattice animal model, Cates and DeutscrP^ 
modified Flory's argument and arrived at the conjecture 
that (R g (N)) ~ TV 2 / 5 for the nonconcatenated unknotted 
rings in the melt. An alternative theory, also based on 
the lattice-of-obstacles picture, but consistent with the 
crumpled globule model for a single ring chairl^, sug- 
gested that rings in the melt are compact objects in the 
sense that (R g (N)) ~ TV 1 / 3 . Lacking any systematic the- 
oretical approaches and based on purely heuristic argu- 
ments, this controversy has remained unresolved. 

Experimental efforts to synthesize melts of pure rings 
have suffered from difficulties in purification and polydis- 
persity. Early studies considered polystyrene and polybu- 
tadiene ring polymers synthesized in both good and theta 
solvent conditions! 41 " 44 ^ However, due to the presence of 
linear chain contaminants and the possible presence of 
knotted conformations, these works have been received 
with skepticism. More recently ring polymers were iso- 
lated from linear chains having roughly the same molecu- 
lar weight with a very high degree of separation ! 45 * 46 ! The 
concentration of linear chains is estimated to be less than 
one chain in a thousand. For these systems, however, 
only dynamical properties have been reported. Kapnistos 
et al. have made stress relaxation measurements for these 
samples as well as samples that were contaminated sys- 
tematically with linear chains! 4 ^ The authors claim that 
the viscoelastic response of the system was significantly 
affected already at concentrations of linear contaminants 
well below the overlap concentration. The effect of linear 
contaminants will be investigated in a separate papeJ^. 
Nam et al. have synthesized and investigated melts of 
low-molecular-weight polydimethylsiloxane rings 4 ^. Syn- 
thesizing and isolating pure, monodisperse, unknotted 
and nonconcatenated rings remains a major challenge 
and thus computer simulations of perfectly controlled 
systems can be very helpful to learn more about these 
melts. 

The early theoretical investigations were followed by 
a series of computational studies EliMtHl i n the first 
worlpo', the bond fluctuation model (BFM) 59 was em- 
ployed, and rings up to N = 512 monomers were exam- 
ined at the density of 0.5 (in units of occupied lattice 
sites). The first results appeared to support the Cates 
and Deutsch conjecture 36 , that (R g ) ~ iV 2 / 5 . However, 
the entanglement length for linear BFM chains at a den- 
sity of 0.5 is known to be fairly large, A^ e? ii near > 70. 
Therefore, the rings examined in this work are still fairly 
short, less than 77V e) ii near , they cannot be expected to 
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be "asymptotic". In the second study from the same 
groupP^the chain length was extended to N = 1024 and 
also stiffer chains were considered, which at the same 
density have a smaller entanglement length. Although 
a crossover to (R g ) ~ TV 1 / 3 was observed, the ultimate 
scaling could not be settled because the overall structure 
of the stiff ring polymers did not resemble a collapsed 
globule, being quite extended to accommodate the higher 
bending energies. More recently some of us studied a 
melt of nonconcatenated flexible rings on a simple cu- 
bic lattice for chain lengths up to N = 1000G 21 . This 
was achieved by a special Monte Carlo method which 
included nonlocal moves to speed-up the relaxation^. 
While the simulations were started from qualitatively dif- 
ferent initial conformations (collapsed globules, extended 
double stranded paths and growing chains) the equilib- 
rium properties were found to be the same. For large 

rings, Vettorel et alPH found the (R^) 1 ^ 2 ~ TV 1 / 3 scaling 
with a wide crossover which can be interpreted as a 2/5 
intermediate "asymptotics" . Suzuki et al.^have also re- 
produced the 1/3 scaling regime. All these works have 
left not only several static properties but essentially all 
dynamic properties unexplored. 

In the first of two papers, we investigate the static 
properties of ring polymer melts using molecular dynam- 
ics (MD) simulation. To our knowledge the only previ- 
ous MD studies of ring polymer melts are by Hur et alP^ 
and Tsolou et alP^who considered polyethylene for short 
chain lengths (i.e., N < 3A/" e) ii near ), and Hur et alP^l for 
chain lengths up to about lliV ej ii near , where AT e5linear is 
the entanglement molecular weight of the linear polymer 
system. Here we consider ring and linear melts composed 
of much longer chains. The number of monomers per ring 
ranges from 100 to 1600 or 3Ar ejlinear to 57N e ^ near The 
linear chain melts were simulated using the same model 
as the rings with chain lengths varying from 100 to 800 
monomers per chain. Such highly entangled melts have 
very long relaxation times and a highly parallelized sim- 
ulation code was used. For the longest rings as many 
as 2048 IBM Blue Gene/P cores were used for a single 
simulation. 



The simulation model is presented in Section II A This 
is followed by a description of the melt preparation and 



the simulation protocol in Section II B The static prop- 
erties for the ring and linear melts are reported in Section 
III followed by a comprehensive discussion of the results 
in Section |IV[ while the key findings of the work and fu- 
ture challenges are presented in the Conclusions section. 
Dynamic properties based on the same simulation runs 
are presented in the subsequent papeffl 



II. SIMULATION METHODOLOGY 
A. Model 

Based on the Kremer-Grest (KG) model 10 , all beads 
interact via a shifted Lennard- Jones potential (WCA or 



Weeks- Chandler- Anderson potential) (Eq. [T]) with a cut- 
off of r c — 2 1 / 6 a. Covalently bonded monomers between 
nearest neighbor beads along the chain also interact via 
the finitely extensible nonlinear elastic (FENE) potential 
(Eq. [2|. Chain stiffness is introduced by an angular po- 
tential (Eq. |3|, where 6i is the angle between bonds con- 
necting beads i — 1, i and z, z + lP^For all ring and linear 
polymer simulations reported in this work: k = 30 e/<r 2 , 
Ro = 1.5 a and fc# = 1.5 e. 
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£4ngle(#i) = k e [l - cos(6»i - 71-)]. 



(3) 



The natural time scale of Eq. Ill is r = cr^/m/e, where m 
is the mass of a monomer. For the simulations we use 
T = 1.0 e/k B and p = 0.85 a -3 . This model has been 
frequently used for simulation studies of polymer melts 
and solutions, so that one can refer to ample information 
throughout the discussion of the present, new results. 
The linear chain data presented in this work is all new as 
well. The ring melts were simulated using 200 polymers 
of length 100 to 1600 monomers per chain. The linear 
systems were composed of 2500 chains of length N = 100, 
250 chains of TV = 200, and 400 chains of length N = 400 
and 800 (cf. Table [i]). An important parameter for our 
discussion is the entanglement length for linear polymers 
in a melt. For the model employed here this parameter 
was determined^ to be 



N e linear — 28 ± 1 



(4) 



Comparing N e ^ neaiT for different simulation models and 
different chemical species allows one to compare equiva- 
lent chain lengths from rather different studies. 



B. Preparation and Production Runs 

The linear chain melts were prepared following the 
procedure outlined in Auhl et al!^ Initially the chains 
are grown as random walks and placed randomly in the 
simulation cell without regard to monomers overlapping. 
The overlaps are removed using a soft potential whose 
strength is increased slowly over a few hundred thousand 
time steps. The Lennard- Jones potential is then turned 
on. The system is further equilibrated using the double- 
bridging algorithm in which new bonds are formed across 
a pair of chains, creating two new chains each substan- 
tially different from the original. This combination of 
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TABLE I. System parameters for the ring and linear melt 
simulations. All simulations were performed with ksT — e 
and pa 3 = 0.85 in a cubic box of side L. The number of 
monomers per chain is N and the number of chains is M. 
The Rouse timeJ^^are r^^mear = N{R^y inear )(/37T 2 kBT and 
Things = tr, iinear/4. The friction coefficient, f, was estimated 
from the initial mean-square displacement data^ of the inner 
monomers of the N — 100 linear system to be 43 r _1 . It is 
found to be 25% lower if the initial ring data are used. 



methods has proven to be a very efficient approach to 
equilibrate long chain melts 

For ring melts the situation is more complex, as we 
do not have well known target conformations. How- 
ever, from the earlier lattice studies we have some insight 
about the general conformations^. The initial configu- 
ration for each ring melt simulation was prepared via the 
following steps^. Rings were constructed with their cen- 
ters at the sites of a simple cubic lattice. The monomers 
of each chain were created on a circular path in a plane 
with a random orientation. By choosing the lattice con- 
stant to be larger than the ring diameter this procedure 
ensures that the rings are unknotted and nonconcate- 
nated. An intrachain attractive interaction was then used 
to collapse the rings. This was done by increasing the 
cutoff distance in Eq. [I]to r c — 2.5 a. The interchain in- 
teractions remain purely repulsive during this phase. For 
the N = 1600 system, an additional force proportional to 
the distance between a monomer and the center-of-mass 
of the chain, was added in order to speed-up the initial 
contraction. 

The density of the system was initially low because 
the rings were widely spaced. A short NPT simulation 
for 500 r was carried out to increase the density of the 
system to the target value (p = 0.85 cr~ 3 ). During this 
phase the applied pressure was slowly increased and the 
intrachain attraction was kept on. Once the target den- 
sity had been reached the NPT ensemble was switched 
to the NVT ensemble. The attractive part of the in- 
trachain bead-bead interaction was switched off and the 
simulations were continued up to about three times the 
Rouse time of each system. The Rouse times for the ring 
systems are given in Table [T[ During this equilibration 
phase the rings diffused a distance of at least their own 
diameter. Taking the typical radius of gyration, in units 
of the entanglement length, the conformations then com- 
pare well to our previous lattice Monte Carlo results on 
very long ring polymers^. In the course of the following 
very long simulations, we also carefully monitored pos- 
sible changes in the instantaneous averages of the span- 



ning diameters and radii of gyration of the rings in order 
to observe any drift pointing towards poor equilibration. 
Within the accuracy of our data, we did not observe any 
drift anymore in all data used for analysis. This is es- 
pecially important for ring polymers starting with some- 
what artificial initial configurations^. 

The production runs were carried out in the 
NVE ense mble with a weak coupling to a Langevin 
thermostat! 10 * 65 * to maintain the temperature. A cubic 
simulation box with periodic boundary conditions in all 
dimensions was used for each simulation. The velocity 
Verlet method was used to numerically integrate New- 
ton's equations of motion. The ring simulations were per- 
formed using ESPResScPI and LAMMP^H with a time 
step and friction coefficient of 0.01 r and 1.0 r _1 , respec- 
tively. LAMMPS 67 was used for the linear systems with 
a time step of 0.012 r and a friction coefficient of 0.5 r _1 . 
In computing mean-square displacements the center-of- 
mass drift of the system was removed. Table [I] gives the 
system parameters for the ring and linear simulations. 



III. CONFORMATIONAL AND STRUCTURAL 
PROPERTIES 

The equilibration of R 2 is shown for the ring systems 
in Fig. [I] The larger chains are found to increase in size 
before fluctuating about an average value while R 2 {t) 
for the smaller chains changes very little over the course 
of the simulation compared to the initial starting con- 
formations. Although long simulation times are needed 
to equilibrate the systems, particularly with large N, our 
runs are long enough because as shown in Fig. [I] even for 
the N = 1600 case, the time window over which average 
quantities were computed is about ten times longer than 
the time required for the system to reach equilibrium. All 
quantities reported in the present and subsequent paper 
are exclusively based on the simulation regime where the 
polymers are fully equilibrated. The time window of this 
regime is approximately five times the longest correla- 
tion time of the N = 1600 ring system, as determined by 
various quantities (cf. subsequent paper on dynamics^). 
The linear systems have longer relaxation times and the 
simulation time of the N = 800 linear system did not 
exceed the longest relaxation time (see below). 

Fig. pj shows the equilibrium average values of R 2 ver- 
sus N for the ring and linear systems. The linear chains 
in the melt exhibit the expected Gaussian behavior or 
(R 2 g ) ~ N. We have to emphasize that this Gaussian be- 
havior is not altered by the recently discovered slowly de- 
caying tangent- tangent correlation along the chain in the 
melt 68 69 . As Fig. [3] demonstrates, the tangent-tangent 
correlation between two points on the chain in the melt 
decays as s~ 3 / 2 with contour distance s, unlike exponen- 
tial decay for the ideal chain in a dilute solutio n. In our 
system, unlike the simulations of Wittmer et al! 68 * 69 *, we 
also observe an initial exponential decay of the correla- 
tion at small s because our chains have some significant 
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FIG. 1. Radius of gyration squared as a function of time for 
the five ring systems. 
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FIG. 2. Equilibrium values of the radius of gyration squared 
(Rg) versus N for the ring and linear systems. The dashed 
lines indicate the crossover and asymptotic scaling regimes of 
the rings. The collapsed ring data was found by simulating 
single rings in vacuum with r c — 2.5 a. The solid line corre- 
sponds to a homogeneous sphere or (R 2 g ) = (3/5)(3N/4np) 2/3 
with pa 3 = 0.85. The error bars for all data points are smaller 
than the symbol sizes. 



bending rigidity (see above). The algebraic decay of the 
tangent-tangent correlation does not alter the conclusion 
that linear chains in the melt are Gaussian in the sense 
of a normal distribution of the end-to-end distances or 
simply the linear scaling of (Rg) with N. This is simply 
because the sum of s -3 / 2 converges at large s, indicating 
that an effective segment comprises only a finite number 
of monomers. 

The ring systems show various regimes depending on 
their lengths. For small N they appear Gaussian. An ex- 
tended crossover regime then follows where (Rg) ~ TV 4 / 5 . 
Finally, the longest chains approach a TV 2 / 3 scaling, show- 
ing that very long ring polymers, with N > 400 ~ 
15A^ e? ii near (see Eq. (Hi)) are required to approach the 



FIG. 3. Tangent-tangent correlation function against contour 
distance between points, s, in double logarithmic scale. The 
correlation function is found to decay much more rapidly for 
the rings (solid line) than the linear chains (dotted line). 
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FIG. 4. Probability distribution function of the radius of gy- 
ration for the (a) rings and (b) linear chains. The insets show 
the probability distribution function of R g normalized by its 
root-mean-square value. 



asymptotic regime. Exactly the same is observed for the 
spanning distance squared R%(N) of the rings, which is 
defined as the mean-square distance between monomers 
1 and N/2 + 1 (or k and k + N/2 with any k) in the 
ring (cf. Table O. These results are in agreement with 
previous simulations^. 

The equilibrium fluctuations of the gyration radius are 
given in Fig. [4] for each of the ring and linear systems. In 
this paragraph we reserve the notation R g for the instan- 
taneous (fluctuating) value of the gyration radius, while 
the mean-square gyration radius is written as (Rg)- The 
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distributions become wider and asymmetrical with in- 
creasing TV, however they follow a perfect scaling when 

normalized by (R 2 ) 1 ^ 2 . The insets to Fig. [4] give the 
probability distribution functions of R g normalized by 
the respective root-mean-square values for the (a) rings 
and (b) linear chains. All five data sets for rings of 
different lengths are found to collapse nearly perfectly 
onto a single curve. For the linear chains, while data 
for N = 100, 200, and 400 collapse about as well as for 
the rings, the distribution for the N = 800 case slightly 
deviates from the others. This is indicative of a system 
that may not have been completely equilibrated despite 
an extremely long simulation time. Some properties are 
more sensitive to improper equilibration than others. For 
this work only the shear relaxation modulus, which is dis- 
cussed in the subsequent paper, is affected by the lack of 
equilibration for the N = 800 linear case. Our data sug- 
gest that all other linear and ring systems are properly 
equilibrated. 

The average values of R 2 and R 2 as well as the ratios 
of the average eigenvalues of the gyration tensor for the 
rings and linear chains are shown in Table [TTJ Not only 
are the rings significantly more compact than the linear 
polymers, they also display a quite different shape. For 
the linear systems the mean-square end-to-end distance is 
found to be related to the mean-square gyration radius 
via (R 2 ) / {R 2 g ) = 6.00 ± 0.02, in excellent agreement 
with the expected value of 6 for ideal chains; this con- 
firms once again that linear chains in the melt are Gaus- 
sian. For the rings, a similar role is played by the ratio of 
the mean-square spanning distance (R 2 ) (mean-square 
value of the vector connecting two monomers N/2 apart 
along the ring) to the mean-square gyration radius (Rg)] 
this ratio for the rings varies from (R 2 ) / (R 2 ) ~ 3.0 for 
N = 100 to about 2.75 for both the N = 800 and 1600, 
indicating that with the approach of (Rg) ~ TV 2 / 3 this 
ratio is close to asymptotic. This is supported by plot- 
ting (Re) I (R 2 g ) versus 1/N and extrapolating to infinite 
chain lengths where one finds a value of 2.75 ± 0.03. For 
Gaussian rings (R 2 e ) / (R 2 ) — 3; while short rings are 
close to Gaussian by this measure, the long ones devi- 
ate noticeably from Gaussian statistics. Both ring and 
linear architectures are found to have shapes that resem- 
ble prolate ellipsoids with the linear chains being signifi- 
cantly more elongated. As TV increases the shape of the 
rings becomes a little closer to spherical with the eigen- 
value ratios of the largest chains being 5.2:2.0:1 (to be 
contrasted to t he eig envalue ratios of roughly 12:3:1 of 
a Gaussian coil) ! 1Q I 7Q I This behavior agrees well with the 
lattice simulations of Vettorel et alP^when the difference 
in A^Knear between the two models is taken into account. 

Further information about rings and chain shapes can 
be obtained from higher moments. Specifically, we can 
define the following series of characteristic lengths: 



In this formula, V{ is the position vector of monomer 
z, while rcM is the position vector of the mass center of 



the entire coil. Obviously, (R 2 g ) 1/2 = Rg Z) (m = 1). Mo- 
ments with higher m are sensitive to extending "penin- 
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sulas" of the structure. Table III shows moments up to 
772 = 5 for both the rings and linear systems. For the 
linear chains, the results are in nearly perfect agreement 
with the theoretical values for Gaussian coils, confirming 
yet again that chains in the melt are practically Gaus- 
sian. By contrast, the situation with the rings is much 
less simple: despite the (approximate) TV 1 / 3 scaling of 
(R g ), characteristic of a compact object, the values of 
the higher moments are nowhere near those for the hard 
sphere, and, surprisingly, not very far from theoretical 
results for Gaussian rings. This indicates that although 
rings in the melt are overall relatively compact, they 
do have many loops which protrude quite far from the 
center-of-mass of the chain. 

To illustrate this Fig. [5] shows three individual ring 
snapshots for each of the N = 100, 400 and 1600 sys- 
tems. Chains are taken from the wings and center of the 
distribution of R g . The extensions increase from left to 
right with the chain in the middle chosen to have approx- 
imately (Rg)- The smallest chains for each system are 
found to be collapsed while the largest are expanded and 
exhibit voids which are occupied by neighboring chains. 
This last point is confirmed by a primitive path analy- 
sis in the subsequent paper^ and self-density data which 
is presented below. The conformations in the left and 
center columns of Fig. [5] are representative of the crum- 
pled globule where each chain has a globule-like core sur- 
rounded by protrusions of size 0(N e ). The variety of con- 
formations in Fig. praises the point of whether there is 
a significant entropy barrier between the quasi-collapsed 
state of a ring and an extended conformation, which re- 
quires significant ring-ring interpenetration without con- 
catenation. To illustrate that this barrier is very small or 
even nonexistent (as suggested by the unimodal distribu- 
tions in Fig. [4^i), Fig. [6] shows R g (t) and snapshots for 
three individual polymers from the N = 1600 ring sim- 
ulation. Chain 35 has a highly extended conformation 
with a large value of R g around t = 34.5 x 10 6 r, how- 
ever, 4.5 x 10 6 r later the chain is found to assume a less 
extended conformation with a corresponding reduction 
in R g of 50%. A similar behavior is found for chain 38 
while the size of chain 51 remains at approximately the 
root-mean-square value (as indicated by the horizontal 
line) over the entire time window. 

The internal structure of the chains may be char- 
acterized by the mean-square internal distances. Let 
[d(s)] 2 denote the average squared distance between two 
monomers separated by s bonds. This quantity nor- 
malized by (R%) is shown in Fig. [T^i. The local rigid- 
ity of the chains leads to an approximate s 2 scaling for 
short separations. An intermediate regime with a lin- 
ear scaling in s is found for lk/d(s = 1) <C s <C N/2, 
where is the Kuhn length which for our modeP^ is 
2.79 a. The [<i(s)] 2 ~ s behavior is found for linear 
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FIG. 5. Individual chain snapshots for three of the ring systems. The left and right columns show small and large chains, 
respectively, while the middle column shows a chain having nearly the average value of R^. 




FIG. 6. Radius of gyration for three individual chains from the N — 1600 ring melt simulation as a function of time along with 
instantaneous configurations. The horizontal line indicates the value of (i?^) 1 ^ 2 . 
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Rings 


Linear 


N 


(1%)/** (Ri)/** (Ai>/(A 3 > (A 2 )/(A 3 ) 


(SZ)/o» (Ri)/* z (Ai)/(A 3 > (A 2 )/(A 3 > 


100 
200 
400 
800 
1600 


17.2 (0.4) 50.8 (1.5) 6.4 2.3 

30.8 (0.7) 88.8 (2.7) 5.9 2.2 

52.9 (1.2) 149.4 (4.8) 5.5 2.1 
87.6 (1.9) 242.4 (7.3) 5.3 2.0 
145.6 (2.8) 401.7 (10.7) 5.2 2.0 


43.4 (1.2) 263.8 (1.6) 12.9 2.8 
88.9 (1.2) 538.9 (1.6) 12.6 2.8 
180.8 (1.3) 1095.3 (1.6) 12.3 2.8 
359.1 (1.3) 2164.9 (1.7) 11.9 2.7 



TABLE II. Size and shape of the ring and linear chains. (R 2 ) is the mean-square gyration radius and (-Re) is the mean-square 
spanning distance between monomers N/2 apart for the rings and the mean-square end-to-end distance for the linear chains. 
Error bars for these quantities are indicated in parentheses. The average eigenvalues of the gyration tensor are arranged as 
(Ai)>(A 2 >>(A 3 >with(^) = E(Ai). 







Linear 


Rings 


Solid 
Sphere 




m 


N = 100 


200 


400 


800 


Gaussian 


N = 100 


200 


400 


800 


1600 


Gaussian 


d(4) /d(2) 

tig /rig 


2 


1.18 


1.18 


1.19 


1.18 


1.19 


1.13 


1.14 


1.14 


1.15 


1.15 


1.14 


1.04 




3 


1.34 


1.35 


1.36 


1.36 


1.36 


1.24 


1.26 


1.27 


1.28 


1.28 


1.25 


1.07 


Itg 1 tig 


4 


1.48 


1.50 


1.51 


1.52 
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1.35 


1.36 


1.38 


1.40 


1.41 


1.36 


1.10 




5 


1.61 


1.64 


1.65 


1.67 


1.67 


1.44 


1.46 


1.49 


1.51 


1.52 


1.46 


1.11 



TABLE III. Average higher moments of the single-chain distance distribution for the linear and ring systems as computed by 
Eq. [5] The values for a Gaussian linear chain and Gaussian ring are also given for comparison with the MD results. For a solid 
sphere R g 2rn) / R g 2) = (5/3) 1/2 [3/(2m + 3)] 1/2m . Note that {R 2 g ) 1/2 = R g 2) . 



chains where it is also the asymptotic scaling. For the 
rings at very large separations a s 2 / 3 regime is found be- 
fore [<i(s)] 2 becomes independent of s as s approaches 
half of the chain. In Fig. [7)3 the mean-square internal 
distances are normalized by those of a Gaussian ring 7 -^ or 
([d(s)} 2 ) = l k l/[l/s ± l/(N - s)], where I is the average 
bond length. Although the overall shape of the d(s) de- 
pendence is quite reminiscent of that of a Gaussian ring, 
it is definitely not the same. Of course, this is not sur- 
prising because unlike linear chains, rings in the melt are 
not Gaussian. 

The static structure function of the individual rings 
multiplied by q 3 is shown in Fig. [sj For 2ir/ (i? 2 ) 1 ^ 2 < 
q <C 2n/lk one expects S(q) ~ q~ x l v for self-similar 
structures. While this scaling is almost perfectly ful- 
filled for linear chains with v — 1/2, the situation is 
more complex for the rings. If the large rings behaved 
like homogeneous, compact objects then one would ex- 
pect S(q) to be dominated by Porod scatterin g 72 * 73 * or 
S(q) ~ q~ 4 . This is certainly not the case here. The 
best fit to S(q) ~ q~ x in this intermediate self-similar 
scaling regime is x = 3.2 ±0.1. This will be explained 
in a more general context in the discussion section. The 

&1 /2 
shows S(q)/N versus q{R 2 g ) • For 

large values of q one finds a gradual increase in S(q)q 3 
indicating a conformation between an overall collapsed, 
but very irregularly and self-similar ly shaped object. All 
systems show a qualitatively self-similar behavior over a 
wide range of q indicating that the general chain shape 
is the same between cases. The shift in the overall am- 
plitude of S(q)q 3 in the intermediate regime, however, 
indicates that, unlike for linear chains, there is some de- 



viation from the ideal fractal object structure, i.e. a weak 
ring length dependence of the density of scatterers. 

Taking the information revealed by Figs. [7|and[8j the 
apparently changing exponents from v = 1/2 for very 
short rings to v = 1/3 for the long ones can be seen as a 
crossover into the asymptotic regime without significant 
qualitative influence on the internal structure. Actually, 
one would expect each ring to have a fairly large number 
of surface beads (n sur f), where a bead is a surface bead 
if it is within a distance r of at least one bead of another 
chain. For random walks, i.e. linear melts, "all" beads 
(O(N)) are surface beads, while for a compact sphere 
only (D(N 2 / 3 ) are surface beads. For the rings we find 
n SU rf ~ P w 0.95 (see Fig. [9]). This power law 
covers the whole range from N = 100 to 1600 without 
any sign of a crossover. The observed exponent of 0.95 
is very close to unity, so one is inclined to suspect that 
in fact n sur f ~ AT 1 , however, the data do suggest that 
this power is slightly below unity and there are also the- 
oretical grounds to expect it to be smaller than 1.0 (see 
Appendix). 

The internal structure of the chains can also be charac- 
terized by P(s), which is the probability of two monomers 
of the same chain being separated by s bonds and a dis- 

P(s) is shown for the 
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tance of less than r. In Fig. 
ring melts for the same two choices of the cutoff distance 
as in Fig. [9j For the larger rings one finds a power law 
of P(s) ~ s~ 7 for intermediate values of s. The value of 
7 is found to decrease weakly with increasing s. For the 
larger cutoff distance and N — 1600, 7 varies from 1.17 
when the data is fit to values of s between 30 — 50 to 1.00 
for values between 110 — 130. The same trend is seen for 
the smaller value of r. 
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FIG. 7. Internal structure of the rings as characterized by the 
mean-square internal distances normalized by (a) the mean- 
square spanning distance and (b) the scaling for a Gaussian 
ring or ([^(s)] 2 ) = W/[l/s + l/ (N— s)], where s is the number 
of bonds separating two monomers and 1 < s < N/2. Nor- 
malizing by (Re) m ( a ) causes all five data sets to coincide 
at s = N/2. 
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To finally confirm that the rings do not assume double- 
stranded conformations we computed the probability dis- 
tribution function of the angle formed by vectors a = 
r^ +n — Vi and b = r J+n — rj, where is the posi- 
tion vector of monomer k and j > i. The calcula- 
tion was performed only when \ri — rj\ was less than 
a distance, r, which was taken as a multiple of the 
excluded- volume interaction cutoff and covered a range 

from nearest-neighbor separations to approximately the 

1/2 

tube diameter (2 (Rg(N e )) ~ 7 a) as based on the lin- 
ear chain system. For all combinations of n = 1, 2, 3 and 
r/r c = 1,2,3,4,5 the distribution function was found to 
be consistent with the vectors having random orienta- 
tions confirming that rings in the melt are not double- 
stranded. Such a finding casts severe doubt on ap- 
proaches which are based on lattice animal arguments. 

The entropy loss associated with the nonconcatenation 
constraint suggests that a ring melt should have a higher 
pressure than a linear melt at the same density and tem- 
perature. This is supported by the finding that a single 
ring in a linear melt (where there is no concatenation 
constraint) has a larger size than in a melt of ringsP^ 
We can estimate the pressure contribution due to the re- 
stricted topology of nonconcatenated rings by following 
the ideas of Ref. 36. According to that work, entropy loss 
due to topological constraints is of the order of the num- 



FIG. 8. The single-chain static structure factor multiplied by 
(a) q 3 for the rings and (b) q 2 (Kratky plot) for the linear 
chains. 
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FIG. 9. Number of surface monomers per chain versus N 
for two choices of the cutoff distance as indicated. The error 
in values of n sur f is approximately 0.5%, significantly smaller 
than the symbol size, suggesting that the observed exponent 
can be distinguished from 1 (dashed line). 
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FIG. 10. P(s) is the probability of two monomers of the same 
chain being separated by s bonds and a distance of less than 
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FIG. 11. Self-excluded monomer density profiles for the rings 
and linear systems. For (a) and (c), r = corresponds to 
the center-of-mass of the central chain while for (b) and (d) 
it corresponds to the center of maximum density. 



ber of rings contacting a given ring, the quantity called 
K\ below. That corresponds to a contribution to the 
pressure which is of the order of ksTKi/Nv, where Nv 
is the volume of a ring. Given that fc#T = e (see Simula- 
tion Methodology section above), and v ^ d 3 , we arrive 
at the estimate of the pressure as K\/N in the units of 
e/cr 3 . Given that K\ is less than 20 (see below), we ar- 
rive for N = 800 at the pressure contribution of less than 
0.025 e/cr 3 . While this estimate is open for criticism, 
since it is based on the method of Ref. 36 which does not 
produce the correct scaling exponent in (R g ) ~ N v the 
deviation is not so large that the total pressure would be 
significantly changed. 

The above estimate of the topological contribution to 
the pressure is compared to the overall pressure which 
was measured by MD simulation and found to be P = 
4.99 ± 0.02 e/cr 3 for N = 800. Thus, the estimated topo- 



FIG. 12. Log-log plot of the self-density as r — >• multiplied 
by N 1 / 2 versus inverse chain length for the rings and linear 
polymers. p S eif around the center-of-mass for linear polymers 
vanishes with iV _1//2 while for the rings it approaches a chain 
length independent constant of p se u(N — >> oo) ~ 0.31 a -3 . 
Note that p se u(r) = p — p ex (r). 



logical contribution is two orders of magnitude smaller 
than the pressure and is roughly the same magnitude as 
the error bars of our measurements. It is therefore not 
surprising that we obtained an indistinguishable pressure 
of 5.00 ± 0.01 e/a 3 for the melt of linear chains at the 
same conditions. Further work is required to detect the 
topological contribution to the pressure and examine its 
exact scaling behavior. 

The overall properties of a single ring in the melt dis- 
play the general scaling of a compact object. However, 
the conformations are open enough in order to allow for 
a significant interpenetration of rings. Such rather ex- 
tended shapes require that the depth of the correlation 
hole, which measures the self-density of the chains, be 
significantly deeper than for linear melts. For linear poly- 
mers the volume of each chain is shared by TV 1 / 2 other 
chains, resulting in a vanishing self-density of the chains, 
i.e. p se if ~ TV -1 / 2 . For a compact object, independent of 
the overall bead density, p se if eventually has to become 
independent of TV. In Fig. 



TTJthe self-excluded density 
profile is shown for the rings and linear chains. This 
profile is computed in the standard way for Figs. [TT^i,c 
except that the origin is taken as the center-of-mass of 
the central chain and the monomers of the central chain 
are ignored in the calculation. For each system the mag- 
nitude of pex in the vicinity of r = does not vanish. This 
arises from the mutual interpenetration of chains and is 
more pronounced for the linear. p ex (r — > 0) is found to 
approach a constant for the longer rings (N > 400). For 
Figs. pT}p,d the center of maximum density of the central 
chain is used. This is found by a binning procedure with 
a bin size of 5 a. The self-density as r —> is shown in 
Fig. [12] where it can be seen that this quantity for the 
rings approaches a constant as N — >• 00. 

The behavior of p e ^{r) is consistent with the ob- 
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FIG. 13. Number of neighboring chains versus N. K\ (r) is the 
average number of chains within a center-of-mass separation 
distance r. 



served decrease in self-density with increasing N, if esti- 

1/2 

mated from (R? g ) • With the self-density being p se u = 

3N/4:7r (R 2 ) 3 ^ 2 (which assumes a homogeneous sphere), 
we find p se if = 0.34 a~ 3 for N = 100 while for the large 
rings (N > 400) this value drops to around 0.23. Since 
the bulk density is 0.85 a -3 this implies that roughly 

1/2 

70% of the monomers within (R 2 g ) of the center-of- 
mass of the long rings come from other chains. This, 
in combination with v = 1/3 means that the number of 
different rings in direct contact with a given ring should 
asymptotically become independent of N. 

The structure of a polymer melt is partly determined 
by the number of neighbors per chain. Let K\(r) de- 
note the average number of chains whose center-of-mass 
is within a distance r of the center-of-mass of a refer- 
ence chain. This quantity is shown in Fig. 13 as a 
function of N for the rings. As chain length increases 

1/2 

K\(r) approaches a constant for both r = (Rg) and 



1 /2 

r = (Riy, which is in agreement with our arguments 
above (i.e., as N — >> oo, v — » 1/3). This is to be con- 

1/2 

trasted with the fact that K\{(R? g ) ) for the melt of 
linear chains grows as iV 1 / 2 . This agrees with the finding 
that the rings in the melt, although collapsed overall, do 
have significant protrusions into one another. The lat- 
tice simulations of Vettorel et alPH gave a limiting value 
of Ki((R g }) ~ 14 as N — » oo which agrees well with 
the value of around 16 for the N — 1600 system in the 
present work. This interplay of (partial) segregation and 
(partial) interpenetration is expected to play an impor- 
tant role for the dynamics of the melts, which is our focus 
in the following paper. 



IV. DISCUSSION 

While the linear chains show a (R 2 g ) ~ N behavior 
for all chain lengths, the rings behave differently. While 
the conjecture of Cates and DeutscfP^ based on a free 
energy argument that the Flory exponent for the radius 
of gyration of a ring polymer melt should be approxi- 
mately 2 / 5 has reportedl y been confirmed by a number of 
simulatiod 2 - 1 ' 50 ' 53 ' 54 ' 58 ' 7 ^ and theoretical 58 studies, it now 
appears that this corresponds to an intermediat e regime 
with the asymptotic Flory exponent being 1/ j 21 * 54 * 5 - 8 -!. 
The present work is in agreement with the scaling laws 
for these two regimes. Vettorel et al! 2 -^ estimate the on- 
set of the asymptotic regime using an argument based on 
Ae, linear- For the scheme where chains form a core which 
squeezes out the other rings, these authors estimate a 
critical ring length N c to be 



N c w p(47r/3) (R 2 e 

.linear (A^linear)) 



3/2 



(6) 



which in our case leads to N c « 2300. Thus it is 
clear that even with our longest rings, which based on 
a comparison of the entanglement length for polystyrene 
would correspond to a molecular mass of more than 
10 6 g/mol, we only would observe the crossover towards 
the (Rg) ~ TV 2 / 3 scaling. While this estimate assumes 
that the rings squeeze each other out more or less com- 
pletely, the data show that this is not the case. The 
incomplete segregation of rings, in which the self-density 
p se if is about 1/3 of the bulk density inside any partic- 
ular coil reduces the expected crossover chain length to 
N c « 700 — 800 in good agreement with our data. This is 
actually a lower limit for 7V C , since the crumpled globules 
are not ideal spherical objects, and it illustrates that one 
needs extremely long rings to unambiguously investigate 
the asymptotic regime. 

Although rings in the melt appear to be compact ob- 
jects in the sense that (Rg) ~ N 2 / 3 at sufficiently large 
A", their compactness is of a non-trivial character. In 
many cases, they are torus (or even pretzel)-shaped; their 
higher moments are much larger than those of a solid 
sphere; and their average self-density, although it does 
not asymptotically depend on A", remains at roughly p/3 
instead of p. These facts indicate that different rings in 
the melt, although reasonably segregated in the scaling 
sense, have nevertheless significant protrusions into one 
another. This can also be seen in the critical exponents 
describing the system. 

In the familiar linear polymer coil, everything is gov- 
erned by the single Flory exponent which describes the 
coil gyration radius. Although the rings in the melt ap- 



pear to exhibit the asymptotic scaling (R 2 g ) 



1/2 



their properties are not completely determined by the 
value of the exponent (y — 1/3), at least not in a triv- 
ial way. There are two important examples of exponents 
that are not directly related to v. The first describes the 
number of "surface monomers" (see above, Section III), 
i.e., monomers of a given polymer that have contacts with 
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other polymers: n sur f ~ we find that our simulations 
yield (3 « 0.95. We argue that the same power (3 should 
describe the number of contacts between two crumples 
of g monomers each, whether belonging to the same ring 
or to two different rings; this number of contacts scales 
as . The second is the 7 exponent which describes the 
looping probability - the probability that two monomers 
of the same ring, a distance s along the chain, will be 
found next to each other in space. The looping probabil- 
ity scales as s~ 7 . Understanding this exponent for rings 
in the melt is of particular interest because this system, 
as mentioned in the introduction section, is expected to 
capture essential features of chromatin packing in chro- 
mosomes. The exponent 7 can be m easured directly for 
chromatin using HiC method j 27 l 3Q t El w e now argue 
that there is a general scaling relation between exponents 
(3 and 7, namely, /3 + 7 = 2 (see Appendix for derivation). 
Our simulation data, presented above (see Figs. 9] and 
10) are closer to (3 + 7 ~ 2.1. It is not completely clear to 
us whether this result agrees with the theoretical predic- 
tion to within the error bars of the simulation data, or 
the deviation is real. Uncertainty in the simulation data 
is thought to arise from finite size effects. Interestingly, 
the simulation result for 7 « 1.1 is also close to the mean 
field prediction 7 = 3^ (see Appendix) as well as to the 
experimental result for chromatin^ 7 « 1. Thus, our 
current understanding is that f3 is slightly smaller than 
unity, while 7 is slightly larger than unity. The former ob- 
servation indicates that there are many contacts between 
crumples, but still not all of the monomers are in con- 
tact, i.e., there is some segregation between crumples; in 
terms of modeling the chromatin, this is consistent with 
the concept of chromosome territories^. The later obser- 
vation is also consistent with the idea that the system is a 
fractal. Indeed, a true fractal on an unlimited spectrum 
of scales cannot be realized with 7 < 1 (because then 
the total number of contacts per monomer, proportional 
to 8~ 7 would diverge). Thus, the observed values 
of (3 and 7 are consistent with the non-trivial segregated 
fractal structure. 



The fact that we observe deviations from a gen- 
erally expected S(q) ~ q~ x l v scaling in the regime 

1/2 

27r/ (R 2 g ) 1 < q < 2n/cr as well as the finding (3 < 1 
can be rationalized by generalizing a recent investigation 
of 2-d polymer melts by Meyer et al. 75 Their analysis 
is based on a rather general theoretical framework of B. 
Duplantiep 7 -^. While the Duplantier theory applies in 2- 
d only, and no 3-d generalization is available (precisely 
because of the difficulties with topological constraints), 
the scaling considerations can be developed for our 3-d 
system in the following way. In a 2-d polymer melt the 
chains form compact objects and partially segregate. The 
situation here is somewhat similar and we can generalize 
this concept to our problem of a melt of nonconcatenated 
rings. In general, with n scatterers (being the whole ring 
or just a part): 



(7) 



where one finds for q — >• 0, S(q) — >• n. For the inter- 
mediate scaling regime, for a fractal structure one ex- 
pects S(q) ~ (q(R g ))~ x . These two arguments imply 
S(q) ~ n(q (R g ))~ x with x to be determined. The total 
scattering intensity I(q) = nS(q) for intermediate val- 
ues of q can only result from the beads at the surface 
of the ring polymer, as these are the only places where 
a scattering probe experiences a contrast. Consequently 
we expect I(q) ~ n^ 3 /q x . Combining all this leads to the 
scaling relation 2 — xv = (3 or x = (2 — (3)/v and a power 
law of 



S(q) 



~ TV 



2-/8)/* 



(8) 



in the intermediate scattering regime. For a random 
walk, where "all" beads are surface beads, (3 = 1 and 
v = 1/2, we recover the well known q~ 2 power law for 
S(q), which also is very well reproduced by our data. For 
the rings the situation is more complex. The best esti- 
mate for (3 is 0.95 as shown in Fig. [9] This would lead 
to a small amplitude shift in S(q) between N = 100 and 
N = 1600 of about 10 — 20%, which indeed is observed. 
In addition, with v = 1/3 we find (2 — (3)/v ~ 3.1 — 3.2 
also in good agreement with the data. Fig. [l4| shows the 
scaling plot for S(q) displaying the overall consistency of 
the data as well as the result that the best fit of (3 for 
the scattering data is 0.93, which once again agrees with 
(3 less than 1. Keeping in mind that very long rings are 
needed to reach the asymptotic regime makes a precise 
estimate of (3 quite difficult. However, the overall con- 
sistency of S(q) and the independent determination of (3 
strongly support the general scheme and (3 < 1. 

For linear chains longer than N e ^ mear , topological con- 
straints manifest themselves in that a "tube defect" (one 
strand leaking out of its tube) becomes entropically un- 
favorable when the involved excursion reaches a length 
larger than L e = b^/N ej n neaiT . In some form, similar be- 
havior is expected in the case of rings as well. After 
the initial Gaussian-like scaling of very short chains, the 
lattice animal structure then emerges as the next natu- 
ral conformation. However, this structure is restricted 
due to the same arguments, which make tube leakage 
very unfavorable. Any extended branch of a lattice an- 
imal would suffer from a severe entropy penalty. This 
restricts the average branch extension to a volume with 
a radius of approximately the tube diameter or to a chain 
length of about 2A r e? ii near . Importantly, effective obsta- 
cles forcing the topological constraints are imposed not 
only by the surrounding chains, but also by the remote 
parts of the same chain. Unlike for a melt of linear 
chains, this is important for rings, because rings in the 
melt are compact and, therefore, the fraction of other 
pieces of the same chain among the spatial neighbors of 
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FIG. 14. The single-chain static structure factor multiplied 
by the general scaling factor g( 2- ^)/^iV 1-/3 with /3 — 0.93 and 
v = 1/3 for the rings with N > 200. 



any monomer is much higher in the ring system than in 
the linear melt. This eventually creates the crumpled 
globule conformation with rather small branches of size 
0(N e ^ linear) ? bulging in and out at any place along the 
contour. A related lattice of obstacles model suggests 
that the nonconcatenated ring in the melt to some extent 
can be thought of as an annealed branched polymer. It 
is annealed in the sense that the material can move from 
one branch to another in the course of thermal motion, 
as Milner and Newhall have recently emphasized 77 . The 
important difference however is that, as the data nicely 
show, the densely-folded regions can open up and allow 
for significant interpenetration. 



V. CONCLUSIONS 



Thus far we have presented a rather detailed study of 
the static properties of a melt of nonconcatenated rings. 
Even though asymptotically the ring extension follows a 
TV 1 / 3 power law, suggesting that the rings have an over- 
all globular conformation, their structure is far from triv- 
ial. Rings still display significant mutual interpenetration 
and the entropic barrier to form fairly large open loops 
is, in fact, so small that they are frequently observed. 
Though overall compact this leads to the fact that almost 
O(N) of the beads are surface beads. These structural 
peculiarities have important consequences for dynamics, 
as will be shown in the subsequent paper 1 . Though 
already quite extensive, the current analysis leaves out 
many properties of relevance to biology, i.e. chromatin 
packing. Quantities specifically motivated by biological 
problems such as inter- and intra-ring contact maps will 
be studied in a separate publication. 



ACKNOWLEDGMENTS 

The authors are grateful to T. Vilgis, T. Vettorel and 
V. Harmandaris for their comments on an early version 
of the manuscript. The ESPResSo development team is 
acknowledged for optimizing the simulation software on 
the IBM Blue Gene/P at the Rechenzentrum Garching 
in Munich, Germany. We thank Donghui Zhang for dis- 
cussions and references relating to experimental studies 
on cyclic polymers. This project was in part funded by 
the Alexander von Humboldt Foundation through a re- 
search grant awarded to AYG. AYG also acknowledges 
the hospitality of the Aspen Center for Physics where 
part of this work was done. WBL acknowledges financial 
support from the Alexander von Humboldt Foundation 
and the Basic Science Research Program through the 
National Research Foundation of Korea (NRF) funded 
by the Ministry of Education, Science and Technology 
(2010-0007886). Additional funding was provided by the 
Multiscale Materials Modeling (MMM) initiative of the 
Max Planck Society. We thank the New Mexico Comput- 
ing Application Center NMCAC for a generous allocation 
of computer time. This work is supported by the Lab- 
oratory Directed Research and Development program at 
Sandia National Laboratories. Sandia is a multiprogram 
laboratory operated by Sandia Corporation, a Lockheed 
Martin Company, for the United States Department of 
Energy under Contract No. DE-AC04-94AL85000. 



Appendix A: Relation between critical exponents of surface 
monomers and loop factor 

The purpose of this appendix is to derive the relation 



between powers /3 and 7, as mentioned in the section IV 
To begin with, recall the definitions of these powers: 

• Exponent 7 describes the loop factor, namely, the 
probability that two monomers a contour distance 
s apart meet each other in space decays with s as 
s~ 



-7 



• Exponent /3 determines the number of "surface 
monomers", n sur f ~ 

To derive the requisite relation between f3 and 7, 
we start by noting that it is equally well applicable to 
the self-similar internal structure of any particular ring 
squeezed between other rings. In this latter case, expo- 
nent /3 describes the number of monomer-monomer con- 
tacts between any two blobs, or crumples. Specifically, 
if two blobs are of some g monomers each, and they are 
in contact (in the sense that the distance between them 
is about their own size), then the number of monomer- 
monomer contacts between them scales as gP . 

With this in mind, consider a ring in terms of N/g 
blobs, g monomers each. Concentrate on pairs of blobs 
which are a distance s monomers, or s/g blobs apart, 
along the chain. There are N/ g such pairs in the chain. 
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Among these pairs some fraction are in contact. This 
fraction scales as a power of contour distance, i.e. it 
scales as (s/g) -7 ; therefore, the total number of blobs in 
contact is (N/g)(s/g)~^ = (Nfs^)g^- 1 . 

Now, let us consider monomer contacts instead of blob 
contacts. First, monomers are not in contact if they be- 
long to non-contacting blobs. Second, given that there 
are g^ monomer contacts per pair of contacting blobs, we 
can find the total number of contacts between monomers: 
(N / s 1 )g 1 ~ 1+ P . Third, we have to realize that what we 
have counted is the number of contacting monomer pairs 
which are a distance s ± g along the chain. The num- 
ber of monomer contacts a distance exactly s apart is at 
least g times smaller, i.e. {N / s 1 )g 1 ~ 2+ ^ . Finally, the 
number of monomer contacts cannot depend on how we 
counted them, that is, it cannot depend on g. Therefore, 
we arrive at 



7 + /3 = 2. 



(Al) 



As an example, consider a Gaussian polymer coil in 
3-d. In this case, 7 = 3/2 and j3 = 1/2, so formula (Al) 
works. In general, for the Gaussian coil in d dimensions, 
we have j3 = 2 — d/2 and 7 = d/2, so again it works. 

As mentioned in the main text, 7 > 1 is required for a 
mathematically rigorous fractal structure, which is self- 
similar over an infinite range of scales. This implies then 
that (3 < 1; in this sense, different crumples of one chain 
or different rings must be segregated from each other, 
meaning that the number of "surface monomers" scales 
at least slower than the total TV. 

Although we have established the general relation be- 
tween exponents j3 and 7, theoretical determination of 
any one of them remains an open challenge. Mean field 
arguments suggest that for any fractal conformation with 
r(s) ~ s v in 3-d, the probability for two monomers to 
meet within a small volume v should go as v/r 3 (s) ~ 
s~ 3l/ , which means 7 = 3v. This gives the familiar result 
£-3/ 2 for the Gaussian coil, while for a crumpled globule 
this yields s -1 , i.e., an impossible 7 = 1. Of course, this 
only means the inapplicability of the mean field argu- 
ment: monomers s apart remain correlated in some yet 
unknown way instead of being independently distributed 
over the volume of the order of r 3 (s). Nevertheless, it 
is worth noting that 7 was found experimentally to be 
very close to unity for human chromosomes 32 , while f3 is 
found very close to unity in the present work. Clearly, 
7 = 1 is impossible for a rigorous mathematical fractal, 
but it might be very close to unity for a real physical sys- 
tem which is only approximately self-similar even if over 
a very wide range of scales. It remains to be understood 
how all these facts are connected to each other. 
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